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On November 16, 1988, the Principal Investigator, Dr. Theodore D. Tarbell, received a 
letter from Dr. Stanley Shawhan of NASA Headquarters stating that the proposal for “An 
Investigation of Solar Active Regions at High Resolution by Balloon Flights of SOUP 
was accepted for a definition phase. The nominal program was to be four years long, 
consisting of definition phase, confirmation for flight, test flight and long duration flight. 
The letter stated very clearly that the Max ’91 Solar Balloon Program was highly cost 
constrained and that the investigation would have to be carried out for the costs proposed 
or be descoped or cancelled. The project was to begin January 1, 1989. 

The flight investigation phase of this contract was cancelled by NASA Headquarters in 
January, 1990, due to budgetary problems, and we were directed to devote our remaining 
resources to a variety of tasks. The remaining tasks, upon which Lockheed, NASA HQ and 
MSFC agreed, were scientific analysis and publication of data from ground-based observing 
with SOUP instrument components, completing the test report on low-voltage piezoelectric 
transducers for use in the image motion compensation system, and submission of a final 
report on the balloon project definition phase. This was done and the contract finally 
expired in April, 1992. 

In January, 1992, four people from LPARL attended the Solar Research Base Enhancement 
Workshop (on other funds). NASA showed renewed interest in balloon flights of optical 
telescopes, and Tarbell and Title contributed to the RBE report chapters on ground- 
and balloon-based observing. The present contract was opened in May, 1992, to extend 
and augment these definition phase activities for another year. This contract continued 
technical studies of the feasability of balloon flights of the the former Spacelab instrument, 
the Solar Optical Universal Polarimeter, with a modem charge-coupled device (CCD) 
camera, to study the structure and evolution of solar active regions at high resolution. In 
particular, different CCD cameras were used at ground-based solar observatories with the 
SOUP filter, to evaluate their performance and collect high resolution images. Some of 
this data has been collected in coordinated observations with the Yohkoh satellite and is 
being analyzed scientifically. 

In the spring of 1992, the SOUP tunable filter instrument was set up in our lab in Palo 
Alto to evaluate the SOHO MDI breadboard CCD camera. This uses a Loral (formerly 
Ford Aerospace) detector built specifically for MDI, a 1024 x 1024 array of 21 micron 
pixels. We used it to take both high resolution and full disk solar images in narrowband 
light at various wavelengths from 5173 to 6768 Angstroms. This detector has insignificant 
fringing in narrowband light, and the gain variations are removed by simple flat-fielding 
using defocussed images. Longitudinal magnetograms and Dopplergrams were made from 
the filtergrams. The images are very pretty cosmetically and have low noise and high 
linearity; the seeing is not very good and so the full disk images are more presentable. 
Since these detectors and cameras are being exhaustively calibrated and space qualified, 
they are good candidates for a future flight instrument. 

In late April, Dick Shine and Darrel Torgerson travelled to La Palma, Canary Islands, 
and set up the SOUP tunable filter with the OSL brassboard CCD camera, which uses 
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a 1024 x 1024 18 micron detector built by Tl-Japan. This camera works well both with 
the tunable filter and with our narrowband K-line filter (3933 Angstroms), showing vastly 
greater blue and near-UV sensitivity than the Loral chip. We observed from early May 
through early July at the Swedish Solar Tower on La Palma, in four shifts: Dick Shine, 
assisted by Igor Zayer of ARC (1 May to 15 May); Ken Topka, assisted by Pete Nisenson 
of Harvard (15 May to 1 June); Zoe Frank, assisted by Louis Strous of Utrecht (1 June 
to 16 June); and Ted Tarbell, assisted by Han Uitenbroek of Harvard, Jose del Toro of 
the Instituto Astrofisica de Canarias, and Goran Scharmer of Stockholm Observatory (16 
June to 6 July). This contract only supported a small fraction of this observing run; most 
of the effort was supported by the OSL contract and by our collaborators. 

At La Palma, our primary goal was to collect high-resolution longitudinal and transverse 
magnetic and H-alpha observations of active regions, simultaneous with X-ray observations 
of the same areas by the Yohkoh satellite. We almost always followed the Yohkoh target 
region and usually stayed in phase with the sunlit parts of its orbit. We obtained not only 
simultaneous observations of flares, but also of quiet sun, flux emergence, sunspots, pre- 
flare buildup, smaller transients (micro-flares, Ellerman bombs, etc.), filament eruptions, 
and pre- and post-flare loops at the limb. Several mornings had very good seeing, and 
many full day movies of the Yohkoh target region were obtained. These contain joint 
observations with the X-ray telescopes of just about every type of dynamic phenomenon 
which occurs in solar active regions, except for a giant flare. In addition, we observed 
during the May 12 rocket flight of the SPDE UV and X-ray instruments built by Marilyn 
Bruner of Lockheed. Noteworthy observations in June and July were M flares on June 8 
and July 3, the last observed in excellent seeing and fully coordinated; excellent sunspot 
sequences on June 16 and July 2; quiet sun and an emerging flux region on June 30; buildup 
and shear of a large active region on June 22-25, leading to an X flare, which unfortunately 
occured just after La Palma and Yohkoh sunset. 

On July 16-17, the instrument was moved from the Swedish Solar Observatory on La 
Palma to the German Vacuum Tower Telescope on the neighboring island of Tenerife. 
Dick Shine, Darrell Torgerson, and Dnyanesh Mathur set up the instrument for solar ob- 
serving from July 19 - August 13. Although some damage to the tunable filter control 
electronics occurred in the transfer, the instrument was been repaired and began observ- 
ing the sun on July 21. Our German collaborators are headed by Wolfgang Schmidt of 
the Kiepenheuer Institute of Freiburg. The observing concentrated on vector magnetic 
measurements and spectroscopy of sunspots, following the Yohkoh target region much of 
the time. A new polarimeter developed by KIS was used with our filter, a prototype 
candidate for the next flight instrument of this type. It should permit much more sen- 
sitive measurements of the transverse magnetic field components, whose measurement is 
considered crucial for estimating energy storage in an active region and forecasting flares. 
In addition, careful measurements of the telescope and instrumental polarization (Mueller 
matrices) were made, and these will be reduced in collaboration with Axel Hoffman, who 
runs the Potsdam vector magnetograph in formerly East Germany. 

The tests in Palo Alto of the Loral and Tl-Japan CCD cameras with the SOUP tunable 
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filter and the successful observing runs at La Palma and Tenerife, in coordination with the 
Yohkoh satellite, produced a very large amount of data, both scientific and calibration. An 
additional CCD test was the use of integrating Cohu cameras with TI sensors to observe 
in the K line (3933 Angstroms) and the G-band (4308 Angstroms). These cameras were 
provided by the Swedish Solar Observatory on La Palma. They can run at normal video 
rates (30 Hz) or integrate for an arbitrary number of video frames, during which time they 
put out a blank video signal. When used with the SSO video digitizer and frame selection 
system, we were able to expose to reasonable light levels with our narrowband tunable K- 
line filter while selecting frames based on the sharpness in the G-band granulation image. 
This image can also be used afterwards to measure the atmospheric distortion in the K- 
line images, for removal by destretching. The sensors and cameras work adequately, with 
reasonable efficiency (typically about 8 frames of integration was used) and cosmetic image 
quality. Thus we collected large amounts of high-resolution K-line imagery simultaneous 
with the usual tunable filter observations. 

Data processing in Palo Alto has proceeded since the observing run, so that nearly all of 
the tunable filter data has undergone the first pass. This includes flat field, dark current, 
and flaw corection, and recording of movies of each wavelength on videodisk. These movies 
are then used for choosing which datasets to study in detail. Some more detailed work 
has been done with data from May 12 (region of the SPDE rocket rocket experiment), 
several filament eruptions of May 13 and June 15-16, and the flares of June 8 and July 
3. These include the 3 July flare movies in AR 7216; sunspot, filament eruption and flare 
observations in AR 7194 on 14-16 June; vector magnetic sets of a sunspot and complex 
neutral line in AR 7186 on 7 June; and movies of remarkable H-alpha loops on the disk 
and above the limb in AR 7205 on 25-26 June. 

Louis Strous, graduate student of Kees Zwaan of Utrecht, spent April through September 
at LPARL. He was reimbursed for living expenses from this contract, in return for his role 
in the observing run and data analysis. He has written a comprehensive set of routines 
for producing magnetic and intensity movies from raw data, with corrections for both 
atmospheric distortion and blurring. The former are removed by destretching, with careful 
separation of solar flows and evolution from seeing-caused distortion. The blurring is 
partially compensated by boosting the high spatial frequencies in each image to match the 
power spectrum of the best image in the series. He applied these techniques to make a 
movie set of a rapidly growing active region. In Utrecht, he is analyzing the flux emergence 
events in these movies and writing his Ph.D. thesis on this data. 

We have begun working with a long-term visitor to our lab, Dr. Jose Carlos del Toro 
Iniesta, of the IAC in Tenerife, on more accurate calibration of longitudinal and vector 
magnetograms. He has radiative transfer codes for computing Stokes profiles of our absorp- 
tion lines using arbitrary model atmospheres. His initial work on Stokes-V diagnostics in 
Fe I 6302.5 using umbral, penumbral, flux tube and quiet photospheric models was shown 
at the AAS meeting. He has continued studying sunspots by making empirical models for 
the temperature and Doppler structure of penumbrae using tunable filter line profiles of 
the non-magnetic line Fe I 5576. 
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Considerable effort has been devoted to sharing data with our other collaborators, training 
them in the idio6yncracies of the instrument and our software. They include H. Uitenbroek, 
B. Noyes, and P. Nisenson of Harvard; E. Zweibel and graduate student B. Hindman of the 
University of Colorado; J. Lawrence of Cal. State Northridge; H. Jones of NSO (Tucson); 
P. Brandt of KIS; R. Rutten of Utrecht; and S. Tsuneta and graduate students T. Shimizu 
and H. Hara of Tokyo. 

Three research papers publishing work which was partly supported by the contract are in 
the Appendices. The first two were begun before this contract but completed during its 
tenure; since they have already appeared in The Astrophysical Journal , only the title and 
abstract are reproduced here. A poster paper was given at the December AGU meeting 
in San Francisco, entitled ’’Simultaneous Observations of Flares in the Solar Photosphere, 
Chromosphere, and Corona from La Palma and Yohkoh,” by Tarbell et al. Three poster 
papers were given at the January A AS Meeting in Phoenix: ’’Comparison of Active Region 
Facular Contrast Measurements to Simple Models,” by Topka and Title; ’’Data Compres- 
sion Experiments with High Resolution Solar Images,” by Shine and Majani (JPL); and 
”From Filtergrams to Physical Atmospheric Magnitudes: A Prospective Diagnostic, by 
del Toro Iniesta, Tarbell and Ruiz Cobo (IAC). Two of these are also included in appen- 
dices. 
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ON THE MAGNETIC AND VELOCITY HELD GEOMETRY OF SIMPLE SUNSPOTS 
Alan M. Title, Zoe A. Frank, Richard A. Shine, Theodore D. Tarbell, and Kenneth P. Topka 

Lockheed Palo Alto Research Laboratory, 3251 Hanover Street, Palo Alto, CA 94304-1191 
GtiRAN SCHARMER 

Swedish Royal Academy of Sciences, Stockholm Observatory, S- 13336 Saitsjobaden, Sweden 

AND 

Wolfgang Schmidt 
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ABSTRACT 

We have observed the disk passage of relatively simple round sunspots using a narrow-band filter and a 
large format CCD detector and have created magnetograms, Dopplergrams, and continuum images nearly 
simultaneously. The spectral resolution of the filter (88,000) allows the construction of “ spectra " for all points 
in the field of view. The changing patterns of the line-of-sight magnetic fields with position around the spots 
and with position of spots on the solar disk show that the mean inclination of the magnetic field increases 
from 45°-50° to 70°-75° across the penumbra and that there is a fluctuation about the mean of about ± 18°. 
These inclination changes occur on the scale of the penumbral fibrils. In addition the direction changes are 
large enough that a substantial fraction of the magnetic field is horizontal (parallel to the solar surface) from 
the mid to outer penumbra. Our Doppler data show that the Evershed flow tends to occur in the regions 
where the magnetic field is horizontal. Taken together these observations suggest that the Evershed flow is 
confined to the regions of horizontal fields, which resolves the old difficulty of nonparallelism of magnetic and 
flow fields in a sunspot. 

We show that a simple sunspot model with azimuthal variations in inclination, but no azimuthal variations 
of field strength, is free from azimuthal Lorentz forces. The meridional currents which arises from the inclina- 
tion variations are parallel to the field lines. This suggests that a cylindrically symmetric magnetostatic 
sunspot model can be perturbed into one with azimuthal variations in inclination with some adjustment in the 
meridional force balance. 

Subject headings: Sun: magnetic fields — sunspots 


1. INTRODUCTION 

Even a simple round isolated sunspot is very complicated structure (Schroter 1971; Gokhale & Zwaan 1972; Cram & Thomas 
1981; Garcia de la Rosa 1987; Kalman 1991). There are intensity inhomogeneities in the umbra and the penumbra, nonsteady 
outflows in the penumbra (Evershed flows), inward motions of bright points along the bright penumbral fibrils, and various 
irregular wave motions in the umbra and penumbra. Given the complexity of the intensity, flow, and wave fields, it would be 
surprising if the magnetic field geometry were not also complex. In this paper we discuss observations of the line-of-sight magnetic 
and velocity field which suggest that even a simple round spot has a mutually interacting spatially complex magnetic and flow 
structure. 

High-resolution magnetic measurements of sunspots are difficult to obtain. Most birefringent filters are not sufficiently narrow to 
resolve spectrally the complex profiles of Zeeman sensitive lines in sunspots. Patrol magnetographs are usually optimized for 
observations outside of spots. High-resolution spectra which measure all four Stokes components with high photometric precision 
often do not have high spatial resolution. Further, vector magnetic field measurements made with filters and spectrographs are 
sensitive to errors introduced by telescope polarization. Even if good data are available, the interpretation is complex because the 
atmosphere of a sunspot is not well understood In particular the local inhomogeneities of the spot may require very detailed 
three-dimensional radiative transfer calculations to interpret spectra properly. 

The above problems have long been recognized. Observers and theoreticians have generated many clever techniques to infer 
properties of spots. These have often used the fact that, as a spot moves across the solar disk, the projection of its vector fields on the 
line of sight changes progressively. In 1908, Evershed (1909) noticed that spot spectra taken with the slit oriented in an east-west 
direction tended to be sloped with respect to lines of constant wavelength. Center (limb) side spectral elements were blue (red) 
shifted. These slopes increased as the spot approached the limb; however, slopes were noticeable even for spots which passed quite 
close to disk center. Evershed inferred that the changing appearance of the spot spectra was due to projection onto the line of sight 
of nearly radial and nearly horizontal outflow in the penumbra of sunspots. Hale et al. (1918) used spectra of Zeeman triplets taken 
in right and left circular polarization to measure where the field was perpendicular to the line of sight, which they called the neutral 
line. By following numerous spots across the disk they were able to generate a mean inclination of the magnetic field as a function of 
radius. 

Hale et al. also observed spots as they crossed the limb in linear polarization at the wavelengths of the n component of Zeeman 
triplets. Because the linear polarization of the n component is in the direction of the field, this technique directly measured the 
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PROPERTIES OF THE SMALLEST SOLAR MAGNETIC ELEMENTS. I. FACULAR CONTRAST 

NEAR SUN CENTER 

K. P. Topka, T. D. Tarbell, and A. M. Title 

Solar and Astrophysics Laboratory, Lockheed Palo Alto Research Laboratory, Dept 9130, B 252, 3251 Hanover Street, Palo Alto, CA 94304 

Received 1991 November 12; accepted 1992 March 9 

ABSTRACT 

We present measurements indicating that the continuum intensity of facular areas in solar active regions, 
outside of sunspots and pores, is less than that of the quiet Sun very near disk center. Analysis shows that the 
observed continuum intensity of faculae at disk center near 5000 A is nearly 3% less than that of the quiet 
Sun The intrinsic value must be even lower than this due to our finite resolution. The continuum contrast 
increases rapidly away from disk center, reaching + 2% at 45°. The zero-crossing point, where the contrast 
changes sign, occurs at 20° heliocentric angle. This is contrary to many earlier observations. The photometric 
data were obtained using the Lockheed tunable filter instrument at the Swedish Solar Observatory, La PalmiL 
Filtergrams of four active regions from 1988, 1989, and 1990 are presented. These were taken with a CCD 
camera and include continuum images of very high spatial resolution (up to 073), and magnetograms (up to 
0745), which have been registered to an accuracy of better than 071. , 

These results are sensitive to the size of the magnetic flux tubes comprising faculae, to their inclinations 
from local vertical at the t( 5000) = 1 level, and to the distribution of sizes of the flux tub cs. The constraint 
these observations place on the size of flux tubes depends upon the value of the zero-crossing point. This is 
important because it is independent of our resolution. Furthermore, because the contrast is changing rapidly 
near disk center, the zero-crossing point is well determined. Our results indicate that most of the flux tubes in 
solar faculae may be very small, in the range 50-100 km in diameter. The results presented here also indicate 
that inclination from local vertical of about 10° at the photosphere is common on the Sun. Footpoints of 
opposite polarity tend to tilt toward one another. 

Subject headings: Sun: facule, flocculi — Sun: magnetic fields 


1. INTRODUCTION 

Within active regions, outside of sunspots and pores, nearly 
all magnetic fields in the low photosphere of the Sun are con- 
centrated into fine structures at or below the resolution limit of 
observations (~T; see Harvey 1977 for a review). Recent 
observations have confirmed and extended these results, 
showing that the smallest magnetic elements are still too small 
to be observed directly, and have field strengths of order 
-2000-2400 G at r = 1 and ~ 1200-1400 G at t = 0.1, with 
diameters of order 100 km (0714, see reviews by Muller 1985 
and Stenflo 1989; see Zirin 1988 for a dissenting view). This has 
led to the concept of the small flux tube, an unresolved region 
of strong magnetic flux surrounded by practically field free 
plasma. Models have been developed for these small flux tubes 
in order to predict their properties (e.g.. Spruit 1976; Deinzer et 
al. 1984a, b; Hasan 1985; Steiner, Pneuman, & Stenflo 1986; 
Knolker, Schiissler, & Weisshaar 1988; Grossmann-Doerth et 
al. 1 988 ; Knolker & Schiissler 1 988). 

Faculae are bright structures seen in the photosphere that 
are associated with magnetic fields. They are especially visible 
in the continuum near the limb, but are also visible at disk 
center on filtergrams taken in photospheric lines. Faculae are 
cospatial with chromospheric plage as seen in the wings of 
strong chromospheric lines such as Ha and Ca ii H and K. At 
very high resolution the faculae consists of many small contin- 
uum bright points, usually smaller than 075, that are also called 
facular points or filigree. Larger structures, including bright 
facular granules and dark magnetic knots, are present in the 
faculae as well. The small bright points can be isolated from 
each other or can be clumped into chains and ringlike struc- 
tures also called “crinkles” (Dunn & Zirker 1973). Isolated 


bright points usually reside in intergranular lanes at the junc- 
tion of three or more granules, and have a lifetime of about 20 
minutes (Muller 1983). All of these observed phenomena, a 
veritable solar zoo, probably owe their existence to the small 
flux tube. Its properties are also important to our understand- 
ing of active region evolution, the interaction between convec- 
tion and magnetic fields, and the heating of the chromosphere 
and corona. 

Title et al. (1990) find, from comparison of very high 
resolution continuum images and accurately registered simul- 
taneous magnetograms, that at the limit of their resolution 
(-0735) and registration (071), the continuum bright points 
(filigree), line center bright points, and magnetic features are 
coincident. The relationship is not one-to-one, however; every 
line center bright point has an associated magnetic feature, but 
many magnetic features do not have associated bright points, 
in either continuum or line center, especially at high magnetic 
filling factors. Title et al. (1992, hereafter Paper I) studied the 
properties of magnetic areas in an active region using simulta- 
neous continuum, line center, Dopplergram, and magneto- 
gram movies over 2 hr long. The observations reported in 
Paper I are very well suited to measure accurately the contin- 
uum intensity contrast of faculae. These data consist of contin- 
uum images (near Ni i 6768 A) of very high resolution of a 
region of dense plage, with simultaneous and accurately regis- 
tered magnetograms. Because these data exist as movies over 2 
hr long, the effects of the 5 minute oscillation on intensity 
measurements can be understood and removed. The results, as 
reported in Paper I, are a contrast of —0.7% at 6 = 14° helio- 
centric angle (ji = cos 8 = 0.970). These results are for one 
image only. 
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ABSTRACT 


We analyse simultaneous, or near-simultaneous, co-registered, digital, photometric 
images of solar photospheric intensity and line-of-sight magnetic field. Images were made 
with the Lockheed tunable filter instrument at the Swedish Solar Observatory, La Palma, 
with the Video Spectra-Spectroheliograph system at the San Fernando Observatory and with 
the new NASA Spectromagnetograph at the National Solar Observatory at Kitt Peak. 

We study the disk-center contrasts of small magnetic elements. While active region 
faculae are dark at disk center quiet sun network features are bright. The field element 
populations of these two kinds of features are quite different. Different contrast center-limb 
functions must be used when estimating their irradiance or luminosity contributions. The 
disk center contrasts of active region faculae are color dependent and indicate a depth effect 
related to the H' opacity of the facular atmosphere. This result is important for calibration 
of monochromatic observations of faculae to bolometric irradiance fluctuations. 

We emphasize the value of cooperative observations among installations whose 
differing strengths are complementary. 
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INTRODUCTION 

The Solar Maximum Mission (SMM) satellite and the onboard Active Cavity 
Radiometer Irradiance Monitor (ACRIM) established that the sun is a variable star. Dips 
in the bolometric irradiance of up to 0.3% matched the rotation of large sunspot groups 
across the solar disk [Willson, et al., 1981]. However, the fluctuations could not be 
attributed entirely to sunspots. The fit was markedly improved by inclusion of the effect of 
active region (AR) faculae [Hudson, 1988; Chapman, et al., 1992]. 

Between 1980 and 1989, the SMM/ACRIM experiment also revealed a few tenths 
of a percent variation of the rotationally averaged irradiance in phase with the 11 year solar 
activity cycle: the sun is slightly brighter on the average when it is more active. This may 
arise from fluctuations of widespread solar magnetic fields not associated with active regions 
[Foukal, et al., 1991]. These fields are found in flux tubes in the quiet sun (QS) "network" 
between the supergranular convection cells over much of the sun’s surface. Like faculae they 
represent enhanced emission, and changes in their numbers could add up to significant 
changes in the solar output. We shall see, however, that AR faculae and QS network fields 
contribute to irradiance fluctuations differently. 

As the sun rotates a feature on its surface is viewed at different angles. The contrast 
of sunspots is not strongly correlated with viewing angle. AR faculae, however, have low 
contrast near disk center but appear brighter near the limb. The nature of this center-limb 
variation is related to the structure of the flux tubes making up the facular elements. For 
example, the model of a facula as an evacuated flux tube with hot walls and a cool floor can 
account for much of the observed contrast variation, especially if allowance is made for the 
effect of a range of flux tube sizes [Spruit, 1976]. 

The center-limb variation of contrast also is used to estimate the contributions of 
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faculae to irradiance changes by means of -prosy" measures like the associated Calcium 
plage area. The fractional change in irradiance can be written [Chapman and Meyer. 1986] 

Al = c B a. \i 4>(»0 c(p) ^ 

Here Cp is a calibration constant. A,, is the observed Ca plage area of a region, 4 >(/*) 
describes the quiet sun photospheric limb darkening, and C(ji) describes the dependence of 
facular contrast on M = cos 0, where 0 is the heliocentric angle between the plage and the 
disk center. The factor n in Eq. (1) accounts for foreshortening as one approaches the limb. 
The facular contribution to solar luminosity fluctuations is the integral of the facular excess 
radiation pattern over the 2 tc steradians of the solar "sky 

Ak = <A)c d a b f 1 4>(|1) c(\i) |1 cfp (2) 

Most empirical models of facular contrasts have been zero or slightly positive at disk 
center [Hoyt and Eddy, 1982; Sofia, et al., 1982; Schatten, et al., 1985; Lawrence and 
Chapman, 1988] although some have used negative disk center contrast [Foukal, 1981]. 
Measurements of facular contrast at the disk center have yielded positive values, though 
these are quite low [Hirayama, et al., 1985; Lawrence 1988] or have been made at 
wavelengths near the H' opacity maximum around 8600 A [Lawrence, et al. 1988, Jones, 
1992]. A recent study [Chapman, et al, 1992] matched photometric measurements of sunspot 
irradiance deficits and of Ca plage areas to ACRIM measurements made while an AR 
crossed the solar disk. This gave a best fit facular contrast function with high contrast near 
the limb but with a disk center contrast of about *1%. This decreases the estimated 

luminosity change of Eq. (2) by about 10%. 

Since 1988 photometric observations have been made with the Lockheed tunable 
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filter instrument at the Swedish Solar Observatory at La Palma. These include very high 
resolution, nearly simultaneous, line-of-sight magnetic field and narrow band, mid-visible, 
intensity images of active regions covering heliocentric angles between 2° (near disk center) 
and 46°. Facular pixels gave contrasts near -3% at disk center, turning positive at a 
heliocentric angle near 20° [Topka, et al., 1992]. New Lockheed observations reported here 
confirm that disk center AR facular contrasts at visible wavelengths are negative. Similar 
San Fernando Observatory (SFO) observations yield the same conclusion. However, 
observations of QS network fields at Lockheed, at SFO, and also at the National Solar 
Observatory at Kitt Peak (NSO/KP) indicate that the associated contrasts are positive at 
disk center. Different contrast functions C(k) apply to the AR and QS cases. 

Observations at Lockheed, at SFO, and at NSO/KP further indicate that disk-center 
facular contrasts are color dependent. This may follow from color dependence of the opacity 
of the solar atmosphere due to the H* ion and may thus be determined by the depth in the 
solar atmosphere to which one is able to see in different colors [Foukal and Duvall, 1985]. 
This is important for modelling facular atmospheres and for calibration of monochromatic 
ground-based observations with bolometric spacecraft observations. 

In the next section we describe the instrumentation and data sets obtained at the 
three observatories which we employ in this paper. Then the results of analysis of the data 
are presented. Finally the results are discussed and conclusions drawn from them. 

OBSERVATIONS 

The Lockheed data were obtained at the Observatorio del Roque de los Muchachos, 
La Palma, Spain, using the 50 cm Vacuum Tower Telescope of the Swedish Solar 
Observatory [Scharmer, et al., 1985] and the Lockheed tunable filter instrument [Title, et 
al., 1992]. This instrument consists of a high-speed steering mirror for image stabilization, 
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reimaging optics, a polarization analyzer, a blocking filter wheel, a narrow bandpass (75 mA) 
tunable filter, and a 1024 x 1024 pixel CCD camera. Only 512 x 512 pixels were recorded, 
however, usually from the center of the CCD. This kept the data cycle to only about 3 sec 
per image. The scale was adjusted in different images from 0.137” /px to 0.166"/px to 
0.274"/px corresponding to fields of view 70", 85” and 140", respectively. Each observing 
sequence had in common the following: one left* and right- circularly polarized pair in 
Fe 1 6302 A for longitudinal magnetograms, a five point line scan (-90, -30, 0, + 30, + 90 mA) 
in Fe I 5576 A or Ni I 6768 A for Dopplergrams, or a longer scan in Fe I 6302 A for the 
determination of Stokes V, and at least one frame in the "continuum." In the Lockheed data 
a search was made within ~ 1 A of the magnetic line center for a spectral interval with full 
width at half maximum (FWHM) = 75 mA as line-free as possible to be used for the 
continuum measurement. For the 6302 Fe line this was at -350 mA from line center - , for the 
green lines it was -750 mA from line center. Exposure times were typically 02 sec. 
Practically, the resolution limit is about twice the pixel scale. For the La Palma observations, 
the resolution for magnetograms was > 0.5" and slightly better for the intensity images. 

The products of interest here are co-registered images of line-of-sight magnetic field 
and of "continuum" intensity. For detailed data reduction procedures see Topka, et al. 
[1992]. Salient features of the data reduction are the following. Time intervals with excellent 
seeing were selected for each observation. The raw CCD images were corrected for pixel-to- 
pixel gain variations, for dark current, and flaws. The magnetograms were formed from the 
blue wing images of the 6302 A line. Because these were not simultaneous with the 
"continuum" images, the two had to be numerically co-registered across the field of view, 
usually to an accuracy of a fraction of a pixel (estimated at 0.1"). The continuum images 
were corrected for any large-scale variations remaining in the flat field. Each image was 
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divided into 1 to 4 sub-images, depending on the nature of the fields in each. The mean 
heliocentric angle is determined for each facular area by comparison to the Mt. Wilson full- 
disk sunspot drawings for that day. The estimated accuracy is 0.25° at disk center and 0.5 
at 40°. Table 1 indicates eight sets of observations included in the present study. 

The SFO data were acquired with the 28 cm Vacuum Telescope and Vacuum 
Spectroheliograph operated in the Video Spectra-Spectroheliograph (VS 2 HG) mode. Details 
of the instrumentation and data reduction may be found in Chapman and Walton, [1989] 
and Lawrence, et al., [1991]. The principal difference between these and the Lockheed data 
is that the CCD images are formed in two circular polarizations with a spectroheliograph, 
so that one CCD axis represents the solar EW direction with scale 0.46 M /px and the other 
axis represents the spectral dimension with scale 8.8 mA/px. The observed spectral range 
included the Fe 1 6302.5 A magnetic line, a telluric 0 2 line at 6302.0 A and a telluric 0 2 line 
at 6302.8 A which is sometimes blended in the umbrae of sunspots. The solar NS direction 
is scanned by moving the spectrograph slit across the solar image. The scale in the NS 
direction is set by the scanning rate and has varied from 0.43"/px to 0.83"/px. The resulting 
data "cube" can be reduced to different kinds of spatial maps by appropriate operations on 
the spectral dimension. A "continuum" image was formed from the average of 5 adjacent 
pixels at 6302.25 A. A line-of-sight magnetic image was formed from the first moment of the 
Stokes V profile of the magnetic line about the center of the I profile; the noise level is 
- 40 G. In this procedure corresponding pixels in the different images are simultaneous and 
co-registered; no de-stretching is required. The effective resolution for these data is > 1". 
The field of view is about 200" EW and variable NS depending on the scan speed and 
duration. Location of features on the solar disk is determined from the suitably rotated 
sunspot drawings in the Solar Geophysical Data Bulletin to an accuracy of 1° near disk 
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center and 2° at 40°. Table 2 indicates seventeen sets of SFO observations included in the 
present work. 

The new NASA/NSO Spectromagnetograph [Jones, et al. 1992] is now in operation 
with the NSO/Kitt Peak Vacuum Telescope [Livingston, et al., 1976]. Data analysis 
techniques associated with it are described by Jones [1992]. The data are similar in concept 
to the SFO data in that a (spatial x spatial x spectral) data cube is formed in two circular 
polarizations, and that these can be reduced to various kinds of spatial maps, such as 
continuum intensity and line-of-sight magnetic field. The ability to choose different spectral 
regions for analysis is much greater than at SFO. Here we have included data in the visible 
at 5507 A and 6122 A and in the near IR at 8688 A. For the Kitt Peak data, the "continuum" 
was taken as the average of the two extreme pixels at each end (i.e. 4 pixels total) of the 
overall spectral interval. The average distance of these from the line center was, for the 
5507 A line, 1.5 times the line FWHM, for the 6122 A line, 1.9 times the line FWHM, and 
for the 8688 A line between 1.4 and 1.9 times the line FWHM. The 8688 A results were 
insensitive to the variation. The pixel scale is 1.14"/px, so the resolution is typically > 2". 
The field of view is the full solar disk. 

NSO/Kitt Peak data included in this study are full-disk continuum and line-of-sight 
magnetic field images made at 6122 A on 22 November 1991; at 8688 A on 7 October 1991, 
16 January 1992 and 3, 4 April 1992; and at 5507 A on 9, 10, 14, 30 June 1992. 

RESULTS 

The contrasts we refer to are the observed percentage excess (or deficit) intensities 
of given image pixels over the average intensity of defined background pixels: 

Q ■ 100 x (Ij- <I>)/<I>. The contrasts are measured in some specified wavelength band 
which we refer to as "continuum." For each of the three data sets discussed, this and the 
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background pixel population are defined differently. Because the chosen wavelength bands 
are not far from magnetic lines, this is not really a continuum contrast. Neither is it white 
light. The goal, of course, is to represent the bolometric irradiance change, so this measure 
should be typical of the spectrum as a whole. A survey of the spectrum indicates a large 

number of lines of various kinds [Kurucz, 1991]. Thus some line contamination of the 

* 

measurement may be to the good (we believe it probably less than 5 %); certainly pure 
continuum would be atypical. Our use of the term "continuum" should be understood with 
this in mind. 

From the magnetic and intensity images, we classify the continuum pixels according 
to their field strength. Field "free" pixels are used to define the QS photospheric level. Then 
intensities are converted to relative contrasts and these plotted versus field strength. Fig. 1 
demonstrates this process for Lockheed observations of three ARs. The average intensity 
in the lowest magnetic bin (*12.5 G to 12.5 G) defines the background level. The facular 
areas at 3° and 17° from disk center show negative contrasts for all field strengths and thus 
have a negative overall contrast. The region at 40° from disk center has positive contrasts 
at most field strengths and overall. 

Fig. 2 depicts a corresponding analysis of SFO data. For illustrative purposes, all data 
points in the range 10° s 0 s 20° were included from the 2nd and 3rd regions in Table 2. 
Shown are contrasts versus line-of-sight magnetic field in 20 G bins, together with the 
fraction of pixels populating each magnetic bin. The pixels are strongly clustered in the first 
three bins (B < 60 G), and we use the weighted average of these bins to determine the QS 
background level. At a magnetic field above 420 G both the pixel contrasts and the magnetic 
bin populations decrease rapidly. For consistent comparison of the Lockheed and SFO 
observations we follow Topka, et al. [1992] in interpreting the pixels beyond this "knee" as 
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representing micropores and in excluding them from calculations of overall contrast. As 
shown in Fig. 1, for the Lockheed data the knee occurs above 600 G. The overall contrast 
for the pixels in Fig. 2 is the weighted average over the range 60G s B s 420 G. This gives 
a slightly positive value in this case. 

Fig. 3 shows the result of several Lockheed determinations of overall facular contrast 
for a number of ARs versus heliocentric angle. Each data point includes magnetic fields in 
the range 188 G < B < 612 G. Below 188 G the points reflect a mixture of B signal with 
non-magnetic noise. At 100 G the mixture is about 50%. Since all the bins in Fig. 1 for 
0 = 3° and 17° have negative contrast, the exclusion of pixels below 188 G cannot change 
the signs of the overall average contrasts. The data points in Fig. 3 represent measurements 
made at different wavelengths, with different image scales and under different conditions 
of seeing and stray light. All of these will affect the magnitudes, but not the signs, of the 
average contrasts to some degree. No corrections have been applied; plots like Figs. 3 and 
4, showing contrast averaged over magnetic field, have this built-in uncertainty. As was 
found in Topka, et al. [1992], the AR contrasts near the disk center are negative, crossing 
to positive at an angle of about 20°. 

Fig. 4 gives a similar analysis of SFO data. The active regions in Table 2 were divided 
into 5° bins of heliocentric angle, and then overall facular contrasts for 60 G < B < 420 G 
were calculated for each image in each bin. Also shown are averages over a consolidation 
of all the images, in 10° bins of 0. A linear regression line has been plotted through these 
points. The dark AR facular contrasts near disk center, found in the Lockheed data, appear 
also in the SFO data. 

We have also observed disk center contrasts associated with scattered QS network 
fields, as opposed to AR faculae. This is illustrated for 0 - 2“ for the Lockheed data in 
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Fig. 5. The contrast versus magnetic field plot starts negative, but then turns positive at 200 
G and remains positive until nearly 500 G. This is markedly different behavior from the 
near-disk-center AR facular contrast curve which is plotted for reference. 

The overall mean contrast is +0.68% for all pixels at 150 G and above. The 
statistical error is ± 0.05%, but the uncertainty due to definition of the background level 
is probably closer to ± 0.2%. This weighted mean is dominated by the peak in contrast 
occurring near 350 G. Above 475 G the contrast is negative due to the presence of 3 
micropores within the field of view. These data, however, total only 59 pixels so they 
contribute little to the weighted average. Between 150 G and 200 G there are nearly 900 
pixels, enough to contribute significantly to the weighted mean, but the contrast is very 
small. 

As shown in Table 2, five observations have been made at SFO of disk center 
network. The contrast versus magnetic field for the aggregate of these observations is shown 
in Fig. 6 along with the magnetic bin occupation fractions. Essentially all of the magnetic 
bins show positive contrast. The overall contrasts for each of the five images are also shown 
in Fig. 4 together with the overall average contrast for all the network pixels. Fig. 6 also 
illustrates the contrasts and bin populations for the aggregate all of the AR observations 
within 10° of the disk center. These give a negative overall contrast represented by the 
leftmost open square point in Fig. 4. The disk center contrast of QS network elements are 
greater than those of AR faculae to a high level of statistical significance. 

Data similar to the Lockheed and the SFO data have been acquired also at NSO/KP 
with the new NASA/NSO Spectromagnetograph. Several sets of full-disk continuum and 
line-of-sight magnetic field maps were made, as listed above, at 8688 A in the near IR. 
These data yielded positive facular contrasts at essentially all magnetic field strengths and 
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at all disk positions. This result matches earlier SFO measurements of positive contrasts at 
8640 and 8680 A [Lawrence, et al., 1988]. NSO/KP observations also have been carried out 
in the mid-visible at 5507 and 6122 A. These give negative contrasts at some magnetic field 
strengths and positive contrasts at others. Fig. 7 presents some of these contrasts for various 
magnetic bins and in 5° bins of 8. In Fig. 8 the data between 0° s 0 s 20° are combined and 
plotted versus line-of-sight magnetic field. These observations parallel the Lockheed and 
SFO observations shown in Figs. 5 and 6. When differences in resolution are taken into 
account, this result is quite consistent with the Lockheed and the SFO observations of disk 
center QS network. The NSO/KP observations were made when no large ARs were near 
the disk center, so they are basically network observations. 

DISCUSSION 

SFO observations of AR faculae yield a contrast versus line-of-sight magnetic field 
dependence which is virtually identical to that of Topka, et al. [1992]. Any discrepancies can 
be attributed to differences in resolution. Negative disk center contrasts were found at SFO 
for facular pixels defined to have field strengths between 60 and 420 Gauss. Analysis of 
Lockheed observations of three additional active regions with heliocentric angles 3°, 17° and 
40° further confirms the dark disk center active region facular contrasts. 

A surprising result comes from observations of disk center, QS network, as opposed 
to AR faculae. One such observation has been carried out with the Lockheed instrument 
and five at SFO. In each case the average magnetic pixel contrast is more positive than the 
corresponding AR contrast. The Lockheed/La Palma case is illustrated in Fig. 5 and a 
composite of the five SFO cases is shown in Fig. 6. Fig. 8 shows a summary of NSO/KP 
observations of the near disk center made at 6122 A on 22 November 1991 and at 5507 A 
on 9, 10, 14 June 1992. These QS network results match those from Lockheed and SFO. 
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These results have a number of implications. First, estimates of irradiance and 
luminosity variability from measured magnetic field areas depend strongly on facular center- 
limb contrast functions. See Eqs. (1) and (2) above. Those who wish to model the facular 
contributions to solar variability both on the active region and the solar qrcle time scales 
[see, for example, Foukal, et al., 1991] may need to use different center-limb contrast 
functions for the AR faculae and QS network. We estimate that differences -100% in the 
estimated irradiance variation could ensue from this. 

Second, Figs. 5 and 6 show that the magnetic pixel distributions for the QS network 
fall off rapidly with magnetic signal and truncate at low field levels: about 500 G for the 
Lockheed data and about 250 G for the SFO data. However, the AR distributions from both 
data sets show long, well-populated tails of dark pixels extending to much higher field 
strengths. These represent the dark features we have called micropores. In the aggregate 
they account for a large fraction of the magnetic flux in the AR faculae. If included in the 
calculation they would make a substantial negative contribution to the overall contrast 

This result exhibits a remarkable difference between the fundamental make-ups of 
the QS network and AR faculae. In particular, micropores are almost completely lacking in 
the network, but are quite numerous in faculae. Micropores are the darkest component of 
AR faculae at disk center, but they may be the brightest component at the limb [Knolker 
and Schiissler, 1988]. Since the network is essentially free of micropores, AR faculae may 
be brighter at the limb, on average, than QS network. 

Our micropores appear identical to the "invisible sunspots" reported by Zirin and 
Wang [1992] and perhaps are related to the "magnetic knots" of Beckers and Schroter 
[1968]. Zirin and Wang reported the presence of invisible spots in network fields far from 
ARs. A few micropores appear in the QS Lockheed image. However, these were localized 
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in a small area of relatively strong field. This may indicate the presence of small, scattered 
patches of "enhanced" network not associated with ARs. 

The brighter disk center emission of QS network versus AR faculae also may stem 

from differing flux tube orientations. The more closely packed AR flux tubes may be more 

closely and more vertically aligned. The more loosely packed QS flux tubes may be observed 

< 

at a variety of inclination angles, and thus of effective center-limb distances. Those seen at 
inclination angles above 20° should appear bright rather than dark. This interpretation is 
supported by the recent observations of Murray [1992], who found inclination angles of up 
to 45° in weakly magnetic areas, but much less in stronger field areas. 

Joint observations like those presented here may reveal different features due to 
differing spatial resolution. For example, the Lockheed observation of the QS network 
(Fig.5) shows a dip in contrast between 0 G and 200 G that is not present in the SFO 
observations (Fig. 6). Statistical analysis indicates that the dip is real, rather than an artifact 
of the Lockheed data. It is present in both the AR and QS network data in Fig. 5. The dip 
must be due to a small and local phenomenon which is resolved in the Lockheed 
observations but not in the SFO observations. This may constrain models of small flux tubes, 
which must be able to correctly predict both the Lockheed and the SFO results. 

A study of facular contrasts near 8688 A [see also Jones 1992] has indicated positive 
contrasts near disk center. This is in agreement with earlier measurements by Lawrence, et 
al. [1988] at 8640 A and 8680 A. These wavelengths are near the H' opacity maximum, and 
their difference from the mid-visible, lower opacity, results may represent a depth effect. 
Facular contrasts near 1.6/* at the H' opacity minimum are dark near disk center [Worden, 
1975; Foukal, et al., 1989, 1990]. This suggests that contrasts at various wavelengths may 
correlate positively with H' opacity. For visible wavelengths, the opacity increases with 
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wavelength, so one can hypothesize that disk center contrast should increase with wavelength 
in the visible range. This is indeed the case for Lockheed observations in three colors of a 
single AR near heliocentric angle 3°. Over a range of magnetic field strengths from 0 to 
600 G, corresponding to facular emission, the contrasts increase in the order of increasing 
wavelength, with the 5576 A contrast nearer the 5250 A than the 6302 A. This is 
demonstrated in Fig. 9 for contrasts corrected for location on the Planck curve. This 
indicates a facular effective temperature about 20 K higher in the red at 6302 A than in the 
green at 5250 A. This in turn indicates that the radiation is emitted from different heights 
in the solar atmosphere. 

The above work indicates the value of measurements at many wavelengths across the 
spectrum. This would allow a bolometric calibration of various monochromatic facular proxy 
measures. At the same time, knowledge of facular contrasts as a function of H* opacity and 
hence atmospheric depth would provide important information for construction of models 


of facular structure. 
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FIGURE CAPTIONS 


1. Average pixel contrast versus line-of-sight magnetogram signal in 25 G bins for all 
image pixels from 3 La Palma observations of active region faculae. The field of view is 
approximately 80" x 80"; the mean heliocentric angle of each observation is 3° (O), 17° (*), 
and 40° (■). Error bars give the 1 o statistical error of each measurement; when not visible 
they are smaller than the symbol. 

2. Average pixel contrast (O) versus line of sight magnetic field in 20 G bins for all 
pixels between heliocentric angles 10° £ 0 s 20° in two SFO images taken at 6302 A on 17 
Aug 1989. Also shown (□) on the same scale is the logarithm of the fraction of pixels in 
each magnetic bin. 

3. Pixel contrast, averaged from 188 G to 612 G, versus heliocentric angle. The smooth 
curve is a quadratic fit to La Palma observations of the 7 active regions listed in Table 1. 
The wavelengths are (O) 5250 and 5576 A, (a) 6768 A and (□) 6302 A. The solid square (■) 
gives the overall contrast of QS network elements near disk center. Statistical error bars are 
smaller than the symbols. 

4. SFO pixel contrast, averaged from 60 G to 420 G, versus heliocentric angle. (O) 
single AR images in 5° bins. (□) averages over all AR images in 10° bins; a linear regression 
line is plotted for these points. (•) QS network contrast measurements. (■) weighted average 
for all network images. Error bars give lo statistical uncertainties when larger than the 
symbols. 
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5. (O) average pixel contrast versus line*of-sight magnetogram signal in 25 G bins for 
faculae in the AR nearest the disk center (3°). (•) contrast versus magnetogram signal in 
50 G bins for disk center network. Logarithms (shifted down by 3 for clarity) of the fractions 
of pixels found in each bin: (□) for AR and (■) for disk center QS network. 

6. (O) average pixel contrast versus line of sight magnetic field in 20 G bins of all SFO 
AR observations with heliocentric angle 0° s 0 s 10° . (•) the same for all five observations 
of QS network near the disk center. Shown on the same scale are the logarithms of the 
fractions of pixels found in each bin for the (□) AR and (■) network cases. 

7. Average pixel contrast versus heliocentric angle in five degree bins for Kitt Peak 
measurements in various magnetic bins: (O) X5507, 4-16 G; (a) A.5507, 16-28 G; (0) A.5507, 
28-100 G; (•) A.6122, 60-160 G; (♦) A.6122, 180-260 G. The 5507 A images were made on 
9, 10, 14 June 1992, and the 6122 A image on 22 November 1991. 

8. Average pixel contrast versus magnetic field in various bins for heliocentric angle 
0 s 20°. (•) data at 6122 A and (O) data at 5507 A. 

9. Average pixel contrast versus line-of-sight magnetogram signal in 25 G bins made at 
(O) 6302 A, (•) 5576 A, and (■) 5250 A. The AR is at 0 = 3°, and the different wavelength 
images are simultaneous within 6 s. If the AR faculae were a perfect blackbody radiator the 
plots would all coincide. 
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TABLE 1. Solar Active Regions Observed at La Palma 


Target . 
NOAA 

6 

Min 

6 6 
Mean Max 

Date 

Time 

(UT> 

Scale 

("/pixel) 

AR 6085 

23 

4.8 

7.4 

06-06-90 

08:25 

0.137 

AR 5168 

10.7 

13.7 

16.7 

09-29-88 

10:07 

0.166 

AR 6086 

21.4 

27.4 

33.7 

06-07-90 

07:57 

0.274 

AR 5572 

37.4 

42.0 

47.0 

07-03-89 

09:15 

0.137 

AR 6774 

0.6 

32 

13 

08-14-91 

10:13 

0.166 

AR 6766 

13.0 

16.7 

20.7 

08-10-91 

11:46 

0.166 

AR 6774 

35.7 

40.0 

44.5 

08-11-91 

10:05 

0.166 

QS 

0.0 

2.0 

3.6 

08-10-91 

17:18 

0.166 
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TABLE 2. Solar Active Regions Observed at SFO 


Target 

NOAA 

6(Deg) 

Min 

6 

Max 

Date 

Time 

(UT) 

Scale("/px) 
NS x EW 


AR 5643 

5 


08-17-89 

18:10 

0.83 x 0.46 


AR 5643 

5 

25 

08-17-89 

18:11 

0.83 x 0.46 


AR 5643 

0 

20 

08-17-89 

22:59 

0.83 x 0.46 


AR 5643 

10 

30 

08-18-89 

22:43 

0.83 x 0.46 


AR 5643 

20 

40 

08-19-89 

20:12 

0.83 x 0.46 


AR 5643 

30 

55 

08-20-89 

11:46 

0.83 x 0.46 


AR 5643 

40 

70 

08-21-89 

20:12 

0.83 x 0.46 


AR 5669 

20 

35 

09-05-89 

19:54 

0.43 x 0.46 


AR 5669 

20 

35 

09-05-89 

19:55 

0.43 x 0.46 


AR 6309 

5 

15 

10-12-90 

20:09 

0.43 x 0.46 


6311 

6312 

AR 6309 

5 

15 

10-12-90 

20:11 

0.43 x 0.46 


6311 

6312 

AR 6734 
6744 

15 

30 

07-26-91 

21:56 

0.43 x 0.46 


QS 

0 

10 

07-06-92 

16:43 

0.43 x 0.46 


QS 

0 

10 

07-06-92 

16:44 

0.43 x 0.46 


QS 

0 

10 

08-22-92 

18:25 

0.43 x 0.46 


QS 

0 

10 

08-22-92 

18:26 

0.43 x 0.46 


QS 

0 

10 

08-22-92 

18:29 

0.43 x 0.46 
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Simultaneous Observations of 
Flares in the Solar Photosphere 
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Palma observing) and NAS8-37334 (SXT). 


1.0 Introduction 


High resolution movies of the photosphere in active repon. reveal fa^inatingflow fields 
and magnetic fine structures which are completely missed m low ^lut.M or mterml 
observations. Simultaneous chromospheric movies m the H«, Na D, and Ca II K 
show a bewildering variety of dynamic features some of which are clearly responses o 
the photospheric drivers. In 1992, we collected very high quality data of h>s type m 
coordinated observations with the coronal X-ray telescopes on Yohkoh. We plan to study 
the best of these movies to identify and characterise quantitatively, when possible, the 
photospheric origins of several well-known beasts in the soo of upper atmospheric activity. 

Excellent observations are available for emerging flux regions, active region growth, flares, 

denies, Ellerman bombs, filaments, and sunspots. For --ft- %L'£S£Z 
tical observation of an M2 flare, with H a , vector magnetic, Na D, and C H K 
of the gradual buildup and impulsive phaseB with 0.5 arcsecond resol . 
simultaneous Yohkoh observations of this entire ewnt in hard and soB :«•»., J^ch 
be coaliined with our movies and compared. Other observations show the buoyant nse 
of tiny magnetic flux ropes through the continuum layer, photosphere low chromosphere 
and eventually (we hope) into the soft x-ray corona. Frequent small flares, with obvious 
reconnections in the H Q fibril topology, are triggered by the emergence and proper mo ion 
of these new flux ropes. 

Our data are time series of very high resolution digital images collected with narrowband 
tunable filters. The 0.05 A SOUP filter isolates several photospheric and chromosphen 
Unefin theTflOO - 7000 A range. A 0.15 A Ca II K fine turmbl. filter 

taneous chromospheric observations. The instrument is operated at two excellent 
telescopes, the Swedish Solar Observatory (SSO) on La Palma and the G "” an V “""” 
Tower Telescope (VTT) on Tenerife. Raw images are converted into coahgned digital 
movies of intensities, Doppler shifts, transverse velocities, and longitudinal magnetograms. 

It is instructive (mid humbling) to recall that raw data not unlike ours, Xtd hi 
movies of many different photospheric and chromosphenc images, were first ccdlecte 
the early 1970’s, with the Universal Birefringent Fdter at Sacramento Peak. , 0f c °"“> 
those movies were on film, and so it proved impo.«ble to compare wavelengths in detail or 
study them quantitatively. The advent of large digital array detectors, fast workstations, 
and Inexpensive data storage has revolutionised high resolution solar physics, 
possible to isolate, identify, and measure physical properties of these complex 
In addition, we also have access to the X-ray movie, of the corona through our Yohkoh 
colleagues. Several of our scientific programs are simply attempts to use these mo ern 
tools and data sets to understand the physics of dynamical features which were discovered 

and named generations ago. 

Since we have more good quality data sets which deserve analysis than people 
to study them, we are eager to discus, collaboration with anyone willing to 

work on the data. 
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2.0 High Resolution Observations 

Starting in August, 1988, we have had an observing run every summer at the 50 ' em 
vacuum refractor of the SSO on La Palma. This telescope is known for its outstanding 
seeing: the highest resolution solar movie ever made is the result of applying LPAKL 
image processing techniques to continuum observations obtained at SSO. Each of these 
five expeditions has collected very high quality time sequences of the type required for 
this proposal. Starting in 1989, we have operated our filtergraph every year in shorter 
runs at the German Vacuum Tower Telescope (VTT) on Tenenfe, while scientists from the 
Kiepenheuer-Insit ut fur Sonnenphysik (KIS) observed simultaneously with their echelle 

spectrograph. 

The purpose of our instrumentation is to collect very high resolution images in narrow 
wavelength bands and well-defined polarization states, with true spatial reso ution down 
to 200 km on the sun. Polarization measurements produce the raw data for longitudin 
and transverse magnetograms. Both vertical and horizontal velocities are measured using 
Dopplergrams and correlation tracking respectively. In the chromosphere, magne ic con- 
nectivity and energy release in flares and other transients are recorded in H a images at 
several wavelength position, within the line. The ground observing instrument consists 
largely of engineering models and flight component, of the SODP instrument and the OSL 

CIP instrument. 

The images are collected with two narrowband tunable filters, one for the mid- visible range 
and a second for the Ca II K line. Each uses a large-format digital CCD camera, with 
images recorded digitally on Exabyte tape. The visible filter system was built initially for 
SOUP and has been upgraded in some components for CIP. The passband is narrow enough 
(50-80 mA) to separate cleanly the core and wings of of Fe I lines m the photosphere for 
magnetic and velocity measurements. The CIP brassboard ^CCD camera ^^th a TI Japan 
uniphase detector is used with this filter. It takes 512x512 or 1024x1024 12 l igi 
images and records them at a rate of one every 2.5 or 8 seconds continuously. 

The Ca II K line filter was built using Lockheed Independent Research funds and is tunable 
across the line with 0.15 or 0.30 A bandpass. It provides chromospheric intensity over a 
range of heights, for studies of wave propagation, chromospheric heating, and mag “e £ 
activity. It can be used with an OSL brassboard camera or with the SSO video CCD 
cameras, which digitize video frames selected for sharpness with optional integration to 
increase light levels. A third camera with a white light or broad band image is also digi lze 
by the SSO system; this is exactly simultaneous with the Ca II image and can be used to 
remove seeing distortions from the narrowband image in subsequent processing. 

The observing sequence is executed by a computer (IBM PC clone), which controls the 
tunable filter system and commands the exposure and readout of the CCD camera. The 
data from the CCD camera is collected on a Vaxstation 3200. The workstation displays 
the filter images on a video monitor and stores the digital images on tape. In 9 , we 

began recording videodisk movies of selected wavelengths in real-time. These movies have 


2 


proven extremely valuable for their instant feedback on evolution of the target region and 
the state of the instrument, as well as for choosing data sets to analyze afterwards. 

A generic observing sequence for active region evolution might take a red continuum im- 
age, a circular polarization pair in Fe I 6302 A for a longitudinal magnetogram, and five 
wavelengths spaced through H a . This could run all day long with a field of view of 140 
arcseconds square, a cadence of about 30 seconds, and a pixel size of 0.27 arcsecond. This 
mode uses the full CCD with 2 by 2 summing. When seeing is especially good, the sum- 
ming can be suspended which yield pixels of 0.13 arcsecond. In order to maintain the 
cadence, the field of view is then usually restricted to the central 512 x 512 area of the 
CCD (field of view 70 arcseconds square). The field of view can be increased to 170 arc 
seconds or 210 arcseconds for very large regions or mediocre seeing. 

In 1992, we spent most of our time at La Palma observing the same region as the Soft 
X-ray Telescope on Yohkoh. For this campaign, we adopted a specific observing strategy. 
During Yohkoh observing times, we ran movie sequences with good time resolution. During 
some of the spacecraft nights, we ran “context” sequences: extensive wavelength scans of 
Ha to show chromospheric connectivity; detailed Stokes polarimetry sequences for vector 
magnetograms; and wavelength scans of circular polarizations for more accurate magnetic 
field strengths. We also attempted to run context sequences in the late afternoon, for 
simultaneous comparisons with magnetograms from US observatories. Different movie 
sequences were available, optimized for different scientific programs such as flux emergence, 
impulsive events, flows at different heights, buildup of magnetic shear, etc. 

The coordinated La Palma- Yohkoh observing was very successful this year. It lasted from 
May 1 to July 6, and several mornings had excellent seeing for periods of an hour or more; 
many full days of movies of the Yohkoh target region were obtained in good seeing. These 
contain joint observations with the x-ray telescopes of just about every type of dynamic 
phenomenon which occurs in solar active regions, except for a giant flare. Noteworthy 
observations were M flares on May 8, May 20, June 8 and July 3, the last observed in 
excellent seeing and fully coordinated; joint observations with the LPARL UV & X-ray 
rocket flight on May 12; excellent sunspot sequences on June 16 and July 2; quiet sun and 
an emerging flux region on June 30; buildup and shear of a large active region on June 
22-25, leading to an X flare, which unfortunately occured just after La Palma and Yohkoh 
sunset. The limitations of solar observing from the ground and low earth orbit are still 
with us! 
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3.0 Flux Emergence and Active Region Growth 

The emergence of new magnetic flux through the photosphere is the first step in the 
formation of an active region. We have several very good observations of young, rapi y 
growing regions. One can learn to recognize emerging flux in movies of raw data and then 
find emergence episodes throughout our movies of mature regions as well. 

We obtained excellent observations on June 30, 1992, of a rapidly growing region near 
disk center. For two hours, the images have very good resolution (one-third arcsecond 
in the best frames). A comprehensive data set was taken, including photosphenc con- 
tinuum, Doppler, and longitudinal magnetic measurements; Na D Doppler and intensity 
measurements; H a line center and wing images; and Ca K filtergrams. In the movies, one 
can see the steady progression of flux tubes rising through the continuum level, middle 
photosphere in the Fe lines, low chromosphere in Na D, and the wings of H a , culminating 
fn the formation of tiny arch filaments in the center of H a . We plan to analyze some of 
these events to make definitive measurements of the buoyant nse of new flux tubes into the 
chromosphere, for comparison with theoretical models of this process, for example those of 
Tajima and his collaborators. Although this was not the Yohkoh target region, sof x-ray 
images are available every ten minutes or so to show the corona at coarser resolution. 
Oddly, this region was not very bright in x-rays in spite of the rapid flux emergence and 

H a activity. 

In every case we have seen, flux emergence does not take place at a neutral line separating 
polarities. Instead, it is occurring in small bipolar events spread over an area large com- 
pared with the size of the bipoles themselves at first appearance. Even a snapshot shows 
that the polarity mixing makes it impossible to draw a neutral line. Perhaps some of the 
rising flux ropes emerge in many places, like multiple loops of a sea serpent above the 
surface; careful study of the magnetic movie may confirm this. Some bipoles also emerge 
with axes at large angles (even right angles) to the main axis of the emergence. The movies 
show that the opposite poles of an emerging bipole separate and then move rapidly across 
the “neutral zone”. Speeds up to 4 km/s are commonly seen in this initial proper motion. 
These exceed the normal horizontal flow speeds measured from granulation, and they ap- 
pear to disagree in direction as weU. These differences indicate the importance of magnetic 
buoyancy and tension forces and they should provide clues to the subsurface configuration 

of the emerging flux rope. 

A very important consequence of this geometry is that many flux cancellation and recon- 
nection events are forced to occur among the newly emerged magnetic elements. With 
opposite polarities moving in opposite directions through the same area, collisions are 
frequently observed. Sometimes dramatic reconnections of the H a fibrils are seen when 
the magnetic footpoints of a set of fibrils are cancelled by intrusion of opposite polarities. 
This is very graphic evidence for magnetic reconnection taking place above the surface, 
presumably in the low chromosphere. These events are discussed further in the section on 

flares below. 
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4.0 Flares 


On July 3, 1992, we obtained superb optical observations in excellent seeing of an M2 
flare in region AR 7216 at N14 E28. The preflare, gradual buildup and 
were fully observed in one Yohkoh orbit, with all instruments performing well and all data 
recovered. Therefore, we will describe these observations in some detail and discuss some 
interesting coronal measurements which can be simply derived from the optical images. 
Figure 1 shows the GOES soft X-ray plot for this event. 

The region was very large (more than 4 arcminutes E-W) and appeared rather simple 
and Inactive in H„. After context observations covering the whole region, we ran a movie 
sequence from 08:55 -09:49 UT; Yohkoh observing started at about 09:00. The movie 
sequence had 3 loops for 3 different time resolutions: H„ hne center and far 'blue ’ wing an 
longitudinal magnetograms every 30 seconds, 5 images spaced through the core of Na D 
evefy 60 seconds, and 5 more H 0 wavelengths and a vector magnetic set every 120 seconds. 
Because of the excellent seeing, the movie used higher magnification, for a 112 
square field (0.22 arcsecond pixels) centered on the neutral line. Figures s 
,2ea in 3 visible wavelengths and soft X-rays at 3 stages in the development of the event. 

At 09:02, a small flare began in a region with emerging mid cancelling flux, in the lower 
left corner of each panel in the figures (Fig. 2). At the same time (comadence or trigger. ), 
weak flare emission began to appear in H„ in the lower right in a purdy monopo ' " ' 

This emission spread in area for the next 40 minute. (Figure 3), following a tortuous 
oath upwards, toward the sunspot in the center. The spreading occured m the following 
P W ay aO flash appeared at the base of each fibril or spicule in the path, followed by 
brightening then disappearance of the individual fibril, dissolving into a diffuse glow in the 
entfre vicinity. This sequence can of course be interpreted in in term, of energy dumped 
suddenly in the low chromosphere at a footpoint, followed by heating mid *^° n .™ k ° h 
did not go into flare mode until 09:37, and the hard x-ray. were not detectable until 09.41. 
At about 09:45, very bright emission appeared m H a and even the wings o a ( PP 
photosphere). At 09:47:30 (Figure 4), the hard x-rays shot up an order of « . 

peaking a minute later (Figure 5). This explosion in hard x-ray. coincide. ^ »»P“* 
jump of optical emission to the pore on the right edge, just above center, of the field. At 
h ”Lm« time there is a sharp increase in H a emission in the far wings, which is s«n n 
very small dot. and arc. on the video images. At 09:49, the La Palma observer .hanged 
back to the large field of view to show the other set of flare footpoints (see F, 8 ur ' ® ’ J 
before Yohkoh sunset at 09:53. After this, the flare cooled gradually and very pretty H a 
loops drained down into both ribbons (Figure 7). 

Neither the H a nor the magnetograms show any apparent cause for this flare: no flux 
collation and no sheared fibrils or neutral line. Both footpoints axe in unipolar regions 
far from any polarity mixing and far from each other. So the instability must have been 
triggered high in the corona. Whatever caused the instability, the gradual buildup to 
it fs^ recorded in the path across the photosphere of emission footpoints. We phmto 
measure the growth of the flare as total footpoint area vs. time throughout the buildup 
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and the impulse phase. Since we have magnetograms accurately realigned with the H 
emission maps, we can convert this curve from area to magnet.e flux as a funct on of 
time The derivative of this, in Mx per second, gives a direct measure of the excitatio 
(reconnection?) rate of field lines in the corona, due to the instability which releases energy 
for the flare. Thus, the high resolution photospheric observations give us an estimate of 
very fundamental coronal quantity. In addition, if we think of the excitation mecharn 
as moving continuously from one field line to another in a limited coronal volume, we ave 
observed the mapping of that motion onto the photosphere. So we have also observed the 
very complicated (perhaps even discontinuous?) path of footpoints of field hues Jong a 
presumably simple continuous path high in the corona. This is relevant to current theones 
of chaotic or turbulent coronal magnetic fields. 

Many small flares in our other movies show clear example, of reconfiguration of the H„ 

fibrils, whose magnetic footpoints are cancelling as opposite polarities collide, i g 

movies, one is compeUed to interpret these as reconnections of the magnetic field hn 

the H 0 chromosphere. The magnetic cancellation may be caused by 

in a neutral sone, as discussed above, or by systematic flow, of swarms of old flux tube 

in a mature region. For example, region, AR 7194 and 7205 each showed ongoing H 0 fibn 

reconnections over periods of several day. in mid-June 1992. Although these flares are not 

impressive in x-ray intensity, they are understandable in term, of their 

and thus they can serve as prototype, for chromospheric reconnection «ents. One could 

select some of the dearest example, of reconnection and to ms*™™* 

average flux cancellation rate, the flux involved in a single H ? fibri event (if this can be 

meaningfully defined), the duration of an event, and speed, of the H Q fibnls. We have no 

idea at This time what the soft X-ray coronal manifestation of these prototype event, will 

be. 
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6.0 Spicules, Ellerman Bombs, and Micro-Flares 


We believe our digital movies are substantially better and certainly more accessible than 
any available in the past. Therefore, there is some hope to make progress in physical 
understanding of these classic problems in chromospheric dynamics. 

Spicules on the disk (or “fine dark mottles”) are seen in the wings of Hq, between about 
0 5 and 1.0 A from line center. They show “noisy” Doppler shifts, with bursts appearing 
and disappearing on timescales of 30 seconds or less. Downflows predominate, but bursts 
of both sign appear in one place with rather different time profiles They seem to nse 
out of the magnetic elements in network boundaries, but the precise relationship with ux 
tubes and granulation remains unclear. We have observed quiet network areas in very 
good seeing with a sequence designed to resolve the motions in spicules. Every 75 secon s, 
it takes photospheric continuum and magnetogram foUowed by 5 sets of 5 H a images each, 
at line center and offsets of ±0.5 and ±0.7 A. Initial viewing of these H« Doppler and 
magnetic movies has suggested the following hypothesis, based on a few clear examples: 
an upward- moving spicule is launched when a small ring of magnetic field (perhaps in the 
lane surrounding a single granule) collapses inwards to a point. We plan to test this more 
systematically by viewing a large sample of events at high magnification, in ^collaboration 
with Y. Suematsu of the National Astronomical Observatory of Japan (NAUJJ. 

The cause of Ellerman bombs is another long-standing mystery. Every active region has 
these emission features seen in the fax wings of H a (not at all in the core), often just outside 
of a sunspot penumbra. They tend to recur in the same place, and their bright ness vaxies 
gradually on timescales of tens of minutes. The movie sequence described above for the July 
3 flare is ideal for studying them, and several others days of high quality are also available. 
Casual viewing of movies and isolated magnetograms reveals small bipoles at the locations 
of some bombs, suggesting the flux cancellation or reconnection models which have een 
published. The proximity to penumbra supports this idea, since nearly horizontal field 
lines are present in the low chromosphere and bends or kinks might reconnect into closed 
loops. Some but not all bombs are also bright in the core of Na D, and strong Doppler 
features may also be present. A competing hypothesis is that bombs are somehow related 
to siphon flows along field fines. We will make sets of vector magnetograms for one or 
more regions with many bombs to see if bipoles are always present. We will also measure 
velocities in the low chromosphere using Na D to search for a hydrodynamical origin. 


At our highest resolution, it seems that every active region always shows tiny, bright, 
explosive events in H a fine center scattered haphazardly throughout. While this may be 
overstated somewhat, transients are extremely common, and a movie viewer often has re 
difficulty deciding whether a region is flaring or not. The exact same comments have een 
made (by Loren Acton, for example) about the SXT movies of active regions, and similar 
transients have been carefully studied in the UVSP data by the MSFC group. We plan to 
coafign a set of H-alpha, magnetic, and soft x-ray movies showing many microflares and 
determine the correspondence (if any) between the chromospheric and coronal transients. 
Are any generalizations possible about the magnetic environment in which they occur. 
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XRAY and OPTICAL DATA BEGIN: 03 JUL 92 
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Figure 1. GOES Soft X-Ray Flux 
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Figure 2. Sample images at the time of the first small flare, around 9:03:40 
on July 3, 1992. Clockwise from top left: Yohkoh/SXT image in soft X-rays, 
longitudinal magnetogram, H a line center, and red continuum. Ground-based 
images have 112 arcsecond square fleld-of-view and 0.5 arcsecond resolution. 
The box on the SXT image shows the La Palma fleld-of-view. 
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Figure 3. Sample images in the gradual phase of the flare, around 9:43:03. One 
of the strands of the twisted filament, which runs vertically on the left side of 
the La Palma field, has just lit up (a micro-flare?). Same format as Fig. 2. 
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Figure 4. The impulsive phase, around 9:48:33, Flare emission has just jumped 
to the pore on the right side. Same format as Fig. 2. 
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Figure 5. Hard X-ray fluxes in 3 energy bands, measure by the Yohkoh HXT 
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Figure 6. Larger fields from La Palma (224 arcseconds square) just after the 
impulsive phase around 9:52:37, showing both flare^ ribbons. 
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Figure 7. Cooling phase around 10:45, showing crossed loops in H a draining 
down into the footpoints. 
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Data Compression Experiments with High Resolution Solar Images 
R. A. Shine (LPARL) and E. E. Majani (JPL) 

Abstract 

High resolution movies are a critical tool for solar physics but they generate enormous 
amounts of data. Data storage and transfers can be a major bottleneck for observations 
and analysis on the ground but it becomes especially acute for solar instruments on small 
spacecraft with limited telemetry and on board storage capabilities. Recently there has 
been a great deal of progress in data compression techniques, especially in the area of 
“lossy” techniques for continuous-tone images. These have the potential of reducing the 
storage and telemetry requirements of images by a factor of 10 or more while maintain- 
ing adequate quality for visual display. We have experimented with lossless and lossy 
techniques for a variety of photospheric, chromospheric, and coronal images. Algorithms 
used include the baseline JPEG (Joint Photographies Expert Group), the EPIC (Efficient 
Pyramid Image Coder), and the HCOMPRESS method. We indeed find that images and 
movies can reproduced from data compressed by factors from 10 to 20 which are virtually 
indistinguishable from the original to the eye. However, a visually acceptable image may 
not be adequate for more detailed analysis of solar movies, such as tracking structures, or 
the combination of images to form magnetograms or dopplergrams. Hence we have also 
experimented with analyzing data with different levels of compression to establish guide- 
lines on the amount of compression tolerable for various types of analysis that might be 
performed on solar movies. 

Methods 

Three publicly available and well documented lossy methods were investigated but after 
preliminary experiments we concentrated on the JPEG algorithm for more detailed studies. 
The 3 techniques are: 

1. baseline JPEG (Joint Photographies Expert Group) algorithm from the Independent 
JPEG Group, archive available on UUNET. This method is based on DCT s of 8x8 
subimages. Compression is based on sorting, quantizing, and run-length/Huffman 
coding of the transform. 

2 . HCOMPRESS, from Richard White at the Space Telescope Science Institute. This 
data compression scheme is being used for the distribution of the digital sky survey 
(done for the guide star catalog for HST) to reduce it to a mere 100 CD-ROMs. It 
uses an H-transform of the entire image. This is a simple 2-D wavelet transform. As 
with the JPEG code, the transform is compressed by quantizing followed by lossless 
encoding of the approximate transform. 

3 EPIC (Efficient Pyramid Image Coder) from the MIT Media Laboratory. This one 
is based on a biorthogonal wavelet (subband) decomposition and a combined run- 
length/Huffman entropy coder. It looks like it makes a tradeoff for a simple decoder 
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at the expense of a more complicated encoder - this is very reasonable for distributing 
images but the opposite of what one wants for spacecraft compression. 

The documentation and sources for these were obtained via anonymous ftp and adapted 
for our image processing environments at both JPL and LPARL. 

All of these methods store a transform of the image and lossy compression is accomplished 
by approximating the transform. Most of the compression comes from approximating 
the higher spatial frequencies (quantizing). They work well (i.e., good reproduction of 
the original with a high compression ratio) when the power is relatively low in the high 
frequency components. This is true of most continuous tone images and many solar images. 

HCOMPRESS can be used in a lossless mode by not quantizing the transform. Similarly, 
JPEG and EPIC can be used in “almost” lossless mode (there is some loss from rounding 
to integers). If the high frequency amplitudes are low, this compression can be significant 
(factor of 4 or more for some solar images). 

The amount of compression (and inversely, image quality) is easily controlled by the level 
of quantization applied. 

Results 

Initial experiments have been done with a variety of non-solar images using all 3 methods 
at JPL by one of us (Majani). Most of the work reported here on solar images has been 

done at LPARL. 

Sample images shown are a x-ray coronal image from the SXT instrument on Yohkoh, a 
photospheric granulation image from Pic du Midi Observatory (digitized from film), and 
filtergrams taken with a CCD at the Swedish Solar Observatory (SSO) on La Palma. We 
have also experimented with other chromospheric and photospheric images (also taken at 
the SSO) with similar results. 


SXT Coronal Image 

Figure 1 shows an image from the SXT instrument on Yohkoh compared with a compressed 
version using JPEG. The original 8 bit image has an entropy of 6.67 bits/pixel and the first 
difference entropy is 4.69 bits/pixel. Our lossless compression code used 4.70 bits/pixel. 
This particular compressed image required an average of 0.75 bits per pixel, more than 
a factor of 10 less than the original. The two are essentially indistinguishable on the 
print and it is even difficult to find differences via blinking them on an image display. 
Table I shows the root mean square (rms) error for a range of JPEG compressions with 
increasing truncation of the DCT amplitudes. To demonstrate how and when the JPEG 
approximation starts to give an unsatisfactory approximation, we have a blowup of a 
portion of the image in figure 2 for a range of approximations. Starting with the x8.4 
compression (corresponding to about 1 bit/pixel) some discrepancies can be seen in some 
of the faint small structures. (These images were printed with a low gamma to enhance the 
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fainter structures where most of the discrepancies are expected.) The problems increase 
with increasing compression but the x20 image is still hard to tell from the original when 
displayed normally. Much more of the high frequency information is lost in the x 53 image 
and the 8 x 8 JPEG cell boundaries are obvious. 

Photospheric Granulation 

Figure 3 shows an image from a series taken at Pic du Midi which was digitized from film 
and then highly processed (this data set is described by Muller et al., 1992, Nature, 356, 
322.). This particular data may be the best granulation movie to date and is representative 
of what could be obtained from a small solar telescope in space. Also in figure 3 are various 
levels of JPEG compression all of which are indistinguishable from the original in this print. 
Figure 4 shows a blowup of a subarea which reveals some differences. Very good quality 
is obtained with the 1.5 bits/pixel image but the 0.75 case is beginning to show some 
distortion of the smallest structures and some of the cell boundaries are beginning to 
become visible. 

Figures 5 and 6 show results using the HCOMPRESS algorithm. It does not perform as well 
as JPEG for this images but this version of HCOMPRESS was “tuned” for photographic 
plates and a version better tuned for the noise characteristics of solar images may work 
better. 

One concern with the JPEG compression scheme was that the 8x8 cells would affect local 
correlation tracking calculations used to study horizontal motions in the solar photosphere. 
For high compression ratios the cell boundaries are readily visible which could cause a sharp 
local peak in the cross correlation at zero lag that might mask or bias the actual offset. 
Numerical experiments with these Pic du Midi granulation images show that any effect is 
small. To further investigate any possible effects, we computed flow maps from a series 
of 100 of these granulation images with 3 levels of compression. The calculations were 
performed in the manner described in the Muller et al. paper. The results were quite good 
and the flow maps from all the compressions displayed in figure 3 gave compatible results. 
Comparing the average flow vectors for the original and the 1.5 bits/pixel compression 
showed an rms error in the velocity amplitudes of about 4% (or about 23 m/s out of the 
564 m/s average velocity) with a bias in the distribution of only 0.4% (or 2 m/s). 


Magnetograms 


The computation of magnetograms and dopplergrams are more sensitive to errors since 
they depend on two or more filtergram images. A special concern for lossy representations 
was that ringing from the truncations of the high frequency components would cause 
artifacts in magnetograms which have regions of high contrast (where the polarity changes) 
and large areas of nearly zero signal. 

Figure 7 shows a magnetogram computed from RCP and LCP filtergrams taken in the 
Fe 6302 line together with the same magnetogram computed from 3 different compressed 
versions of the data. These were done with a 12 bit JPEG code rather than the 8 bit version 


3 


used for the earlier examples. Each was constructed via our usual method of destretching 
the LCP image relative to the RCP image and then computing a line of sight magnetic 
signal from the difference of each pixel pair divided by the sum. The destretching was done 
after the compression and represents another step where compression errors could creep 
in. 


From this example it seems that the filtergrams used to construct magnetograms can be 
compressed down to the 1.0 bit/pixel level without generating any distracting artifacts. 
Table II summarizes some of the errors in the compression of the original filtergrams 
and the rms errors in gauss for the (roughly) calibrated field strength. As would be 
expected, the largest errors are in the dark sunspot umbra where this method of measuring 
magnetic fields breaks down anyhow. The quality is surprisingly good and compressed 
filtergrams can certainly be used to construct magnetograms for most purposes. Some 
types of analysis, such as averaging a large number of magnetograms to improve signal to 
noise in a search for weak structures, are probably not feasible with compressed data. 

Conclusions 

Data compression techniques are now available that can provide an order of magnitude 
reduction in the storage and transmission requirements of solar images with sufficient accu- 
racy for most scientific purposes. This study is not exhaustive but it shows that the JPEG 
scheme is already a feasible technique for spacecraft telemetry and data storage of coronal, 
chromospheric, and photospheric images. Improved data compression methods specialized 
for a particular class of solar images can probably be developed that perform somewhat 
better. The amount of compression can be varied to allow a tradeoff between compression 
and quality for different applications. In general we find that about 1 bit/pixel is sufficient 
to maintain data quality for most applications. Some coronal and chromospheric images 
can be compressed somewhat more (down to the 0.5 - 0.75 bit/pixel range) while well 
exposed photospheric images that will be used for constructing magnetograms or doppler- 
grams may need 1.5 bits/pixel. Effects from the JPEG cell boundaries and ringing from 
the high frequency approximations are surprisingly low for these compression levels. Some 
information is, of course, inevitably loss. In particular, accurate measurements of small 
faint structures will suffer. Also, it is probably not feasible to sum up a large number of 
compressed images to increase sensitivity - e.g., to look for very weak magnetic structures. 

With the current emphasis on smaller scientific satellites with severe telemetry bottlenecks 
(e.g., a Gigabit or two per day for a SMEX), the feasibility of xlO data compression makes 
an enormous difference for solar imaging experiments. 
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Figure 1. An SXT image from Yohkoh taken on 12-May-1992. The upper image is the 
original 8 bit image and the lower was reconstituted from a JPEG compressed file that 
used an average of 0.75 bits per pixel. 


Figure 2. A portion of the SXT image in figure 1 showing the original and 5 levels of 
compression using JPEG. 


Figure 3. A photospheric image of granulation and the image reconstructed from 3 different 
levels of JPEG compression. The bits/pixel required for each is indicated. 


Figure 4. A blowup of part of the image shown in figure 3. 

Figure 5. Same as figure 3 except that the HCOMPRESS algorithm was used. The 
bits/pixel required for each of the 3 compression levels is indicated. 

Figure 6. A blowup of part of the image shown in figure 5. 

Figure 7. Magnetograms constructed from right and left circularly polarized filtergrams 
in the Fe 6302 line for the original data and 3 different JPEG compressions as indicated. 
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Table I. Compression Statistics for SXT Image 


bits/pixel 

compression factor 

rms error (counts) 

rms/mean 

4.59 

1.7 

0.57 

0.008 

2.84 

2.8 

1.91 

0.028 

1.73 

4.6 

2.89 

0.042 

1.21 

6.6 

3.39 

0.049 

0.95 

8.4 

3.68 

0.054 

0.77 

10. 

3.86 

0.056 

0.66 

12. 

3.99 

0.058 

0.50 

16. 

4.30 

0.061 

0.40 

20. 

4.36 

0.064 

0.31 

26. 

4.54 

0.066 

0.15 

53. 

5.16 

0.075 

0.08 

95. 

6.27 

0.091 
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Table II. Compression Statistics for Magnetogram 


bits/pixel 

RCP 

rms/mean 

LCP 

rms/mean 

magnetogram 
rms (gauss) 

2.92 

.002 

.002 

13 

2.06 

.003 

.003 

18 

1.25 

.004 

.004 

28 

0.85 

.005 

.005 

35 

0.52 

.007 

.007 

48 























Comparison of subareas from an SXT image with various JPEG compression factors. These 
are 100x100 pixel subareas of a region showing more structure than typical. The factors apply 
to the entire image. 
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Figure 3. A photospheric image of granulation and the image reconstructed from 3 different levels 
of JPEG compression. The bits/pixel required for each is indicated. 
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Figure 4. A blowup of part of the image shown in figure 3. 
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Figure 5. Same as figure 3 except that the HCOMPRESS algorithm was used. The bits/pixel 
required for each of the 3 compression levels is indicated. 
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Figure 6. A blowup of part of the image shown in figure 5. 
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gure 7. Magnetograms constructed from right and left circularly polarized filtergrams in the 
: C302 line for the original data and 3 different JPEG compressions as indicated. 
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